CN109783972A - Monitoring method based on the internal drainage of check valve leakage current amount that wind-structure interaction calculates - Google Patents
Monitoring method based on the internal drainage of check valve leakage current amount that wind-structure interaction calculates Download PDFInfo
- Publication number
- CN109783972A CN109783972A CN201910089743.3A CN201910089743A CN109783972A CN 109783972 A CN109783972 A CN 109783972A CN 201910089743 A CN201910089743 A CN 201910089743A CN 109783972 A CN109783972 A CN 109783972A
- Authority
- CN
- China
- Prior art keywords
- temperature
- valve
- leak
- hot
- leakage
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Examining Or Testing Airtightness (AREA)
Abstract
The present invention is the monitoring method of the internal drainage of check valve leakage current amount calculated based on wind-structure interaction, the geometrical model that this method establishes pipeline to be analyzed directly carries out fluid and structural simulation according to flow boundary condition, pass through the calculating to branch pipe leak fluid temperature case after supervisor's channel temp, flow and valve before different valves, the outside wall temperature under different internal leakage flows is acquired, valve leakage outside wall temperature characteristic parameter is established, is responsible for valve leakage quantitative prediction model complicated between channel temp and flow, leak fluid temperature and internal leakage flow.Thermocouple is set at the temperature monitoring position of setting, according to the temperature monitoring data temperature characteristic parameter of different circumferential positions, the operating parameter that can be provided according to power plant obtains internal drainage of check valve leakage current amount.This method can be used for core pipeline peace note pipeline internal drainage of check valve leakage current amount monitoring, solve check-valves internal leakage flow information in prior art center pipeline can not quantify and can not accurate evaluation the problems such as.
Description
Technical field
The present invention relates to pipe line check valve internal leakages to monitor field, refers in particular to a kind of based on fluid structurecoupling/temperature field analysis
The monitoring method of internal drainage of check valve leakage current amount is especially adapted for use in the monitoring to core pipeline pipeline internal drainage of check valve leakage current amount.
Background technique
In recent years, in some nuclear power stations, the pipeline failure due to caused by thermally stratified layer is had been reported.A kind of typical heat point
Layer is caused by leaking from RCS (reactor coolant loop) to the turbulent flow infiltration of peace note pipeline and branch pipe internal drainage of check valve.It lets out
Accurately determining for leakage current amount is most important for the Thermal Fatigue due to peace note pipeline, and the effective leakage prison of current still shortage
Survey means.
Method about the on-line checking of industry spot valve leakage at present, mainly has: pressure decline method, vibratory drilling method harmony hair
Penetrate detection method etc..Pressure decline method by respectively installing a pressure sensor in valve upstream and downstream, if upstream pressure be substantially reduced and under
Pressure increase is swum, then can determine that the leakproofness of valve is not good enough.The disadvantage is that needing to punch installation pressure sensor, it is not a kind of
Proper lossless detection method hardly results in popularization and application in practical applications.Vibratory drilling method is to utilize acceleration sensing
Device detects caused pipe vibration signal when valve leakage, the disadvantage is that it is to valve leakage and ambient noise with same quick
Sensitivity, anti-interference ability is not strong, and in the case where tapping pipeline, system is also judged as valve leakage, have it is very high not
Certainty.Acoustic emission detection method is the time and frequency domain characteristics of the acoustic emission signal generated by valve leakage, and detection process compares
Fast, and what does not have on the integrality of valve influences;Disadvantage is that, when internal leakage occurs for valve, the sound that inside generates is sent out
Penetrate that signal component is very complicated, it is difficult for effective identification of the acoustic emission signal of internal leakage, and acoustic emission signal
Sample frequency is higher, this characteristic does not require nothing more than hardware device with very high processing speed and biggish memory space, and
Also there is large effect for the analysis processing of signal, feature extraction and diagnosis modeling.In Chinese Patent Application No.
Disclose a kind of therrmodynamic system valve leakage flow monitoring method in 201611019372.4, the monitoring method include: by
Therrmodynamic system pipeline is divided into multiple infinitesimal pipeline sections, carries out radial and axial build to the infinitesimal pipeline section according to heat transfer principle
Mould, and solved using iterative method along Working fluid flow direction;Set up the large sample model between temperature field and characteristic parameter.However,
More demanding for the pipe of system, the arrangement of valve in such detection method based on temperature, applicability is poor.In short, existing
Perhaps monitoring effect is bad or system structure complexity, detection device are invested greatly for valve internal leak test method, is unsuitable for
More to quantity, spatial distribution complexity core pipeline valve detects.
Since core pipeline pacifies the particularity of note pipeline fluid flow, the prior art uses only the heat transfer of leak fluid
Determine single-point tube wall temperature correlation of leakage flow, this far can not be solved, and there are cold fluid and hot fluid thermally stratified layer, turbulent flows to penetrate into, cold and hot friendship
The problem of for equal complex fluids flowing and phenomenon of Heat, for the disadvantage in above-mentioned art methods, the present invention proposes one kind
Based on the monitoring method of fluid structurecoupling/temperature field analysis core pipeline pipeline valve internal leakage flow, fluid flowing is comprehensively considered
The relevant operating parameter of state determines monitoring position pipeline outer wall temperature leakage characteristic parameters based on fluid and structural simulation, can be real
The quantitative on-line checking of existing core pipe line check valve leakage flow is timely discovery and the maintenance provider of core pipe leakage accident
Just, equally have great importance to the fatigue evaluation of pipeline, and monitoring method is easy to use, engineering calculation requirement can be met.
Summary of the invention
Present invention aims at the demands for core pipeline peace note pipeline internal drainage of check valve leakage current amount Quantitative Monitoring, propose
A kind of monitoring method of the internal drainage of check valve leakage current amount calculated based on wind-structure interaction.This method can be used for core pipeline peace note pipe
The monitoring of line internal drainage of check valve leakage current amount, solving check-valves internal leakage flow information in prior art center pipeline can not measure
Change, generate in pipeline heat fatigue caused by thermally stratified layer, hot and cold alternation etc. can not accurate evaluation the problems such as.
To achieve the above object, the invention provides the following technical scheme:
A kind of monitoring method of the internal drainage of check valve leakage current amount calculated based on wind-structure interaction, is included the following steps:
The first step, the geometrical model for establishing pipeline to be analyzed, and mesh generation is carried out to fluid domain and solid domain, with meter
Fluid operator machine software carries out fluid and structural simulation, obtains valve preceding pipeline temperature field:
Channel temp t is responsible for before known pipeline position to be analyzed valvehot, main pipeline flow Q before valvehotAnd pipe internal leakage stream after valve
The temperature t of bodycoldOpereating specificationDifferent leakage flow QLeak
Tube wall temperature data before lower valve, the corresponding one group of tube wall temperature of each leakage flow;
Based on fluid and structural simulation, monitoring cross section temperature measuring point before n valve in the circumferential direction, the 1st temperature measuring point are determined before valve
Set on the top side of sectional position, the n-th temperature measuring point is located at the least significant end of sectional position;The temperature of two neighboring temperature measuring point successively takes
Average value is denoted as Ti, i=1,2,3 ..., n-1, the area of monitoring cross section is denoted as S between two adjacent temperature measuring pointsi, calculated according to formula (8)
Obtain the equivalent face mean temperature T of tested pointA, maximum temperature T in whole cross section is calculated according to formula (9)topWith minimum temperature
TBottomTemperature difference Δ T, the equivalent face mean temperature T of tested pointAWith temperature difference Δ T as tube wall temperature characteristic parameter,
Δ T=Ttop-TBottom (9)
In formula:
The determination of fluid and structural simulation and leakage flow correlation under second step, different flox conditions:
To fluid and structural simulation is carried out under different flox conditions, tube wall temperature feature before valve is joined using least square method
Number data are analyzed and are handled, and tube wall temperature characteristic parameter is with leakage flow Q before obtaining valveLeak, leak fluid temperature
tcold, supervisor channel temp thotAnd main pipeline flow QhotChanging rule;It fits before valve tube wall temperature characteristic parameter and lets out
Leakage current amount QLeakQuantitative relation formula be T=f (thot,Qhot,tcold,QLeak), T represents T in formulaA,ΔT;
According to T obtained above and QLeakMathematical relationship, solve leakage flow QLeakJoin with tube wall temperature feature before valve
Number, supervisor channel temp thotAnd main pipeline flow Qhot, leak fluid temperature tcoldBetween relational expression be QLeak=g (thot,Qhot,
tcold, T),
QLeak=g (thot,Qhot,tcold, T) and it is tube wall temperature characteristic parameter, supervisor channel temp t before valvehotAnd supervisor
Road flow Qhot, leak fluid temperature tcoldWith leakage flow QLeakBetween complicated valve leakage quantitative prediction model;
Third step carries out real-time monitoring according to prediction model:
Temperature measuring point installs thermocouple before valve, and temperature at temperature measuring point before valve is sent to monitoring server end in real time, is obtained
Thermally stratified layer temperature profile parameter, calculate leakage flow further according to valve leakage quantitative prediction model, compare by two temperature
Both the correspondence leakage flow that degree characteristic parameter obtains, if the leakage flow relative error of the two is less than leak threshold, take
Arithmetic average, if being not less than leak threshold, taking the two the larger value is final leakage stream magnitude, realizes real-time monitoring.
The monitoring method of the above-mentioned internal drainage of check valve leakage current amount calculated based on wind-structure interaction is applied to core pipeline peace note
The monitoring of pipeline internal drainage of check valve leakage current amount.
One or more technical solution provided by the present invention, has at least the following technical effects or advantages:
The invention proposes a kind of monitoring methods of internal drainage of check valve leakage current amount calculated based on wind-structure interaction.The party
The geometrical model that method establishes pipeline to be analyzed directly carries out fluid and structural simulation according to flow boundary condition, before to different valves
The calculating of branch pipe leak fluid temperature case, acquires the outer wall temperature under different internal leakage flows after supervisor's channel temp, flow and valve
Degree, establish valve leakage outside wall temperature characteristic parameter, supervisor channel temp and flow, leak fluid temperature and internal leakage flow it
Between complicated valve leakage quantitative prediction model.It, can be according to the arrangement of pipeline and flowing side in actual field measurement process
Boundary's condition determines that internal drainage of check valve leaks flux prediction model, thermocouple is arranged at the temperature monitoring position of setting, according to different circles
The temperature monitoring data temperature characteristic parameter of all positions, operating parameter (the main pipeline fluid stream that can be provided according to power plant
Amount and temperature, leak fluid temperature) obtain internal drainage of check valve leakage current amount.Error analysis shows this method internal leakage volume forecasting
Precision is high, and method is simple and feasible, can satisfy the demand of the calculating of peace note pipeline heat fatigue, internal drainage of check valve leakage Monitoring and maintenance.
The present invention quantitatively determines internal drainage of check valve leakage current amount by temperature field using the method for fluid structurecoupling, and applicability is more
Extensively.
Detailed description of the invention
Purpose or technical solution in the prior art in order to illustrate the embodiments of the present invention more clearly, below will be to implementation
Example or attached drawing needed to be used in the description of the prior art are briefly described, it is clear that described embodiment is the present invention one
Section Example, instead of all the embodiments.The accompanying drawings in the following description is only the embodiment of the present invention, for this field
For those of ordinary skill, do not make make the creative labor under the premise of, it can also be obtained according to the attached drawing of offer
His attached drawing.In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, and those of ordinary skill in the art are not having
Every other embodiment obtained under the premise of creative work is made, shall fall within the protection scope of the present invention.
Examples of the present invention will be described by way of reference to the accompanying drawings.
Fig. 1 is that check-valves leaks heat fatigue region in the application case study on implementation.
Fig. 2 is heat peace note pipeline geometrical model in the application case study on implementation.
Fig. 3 is grid dividing schematic diagram in the application case study on implementation.
Fig. 4 (a) is the Temperature Distribution cloud atlas of overall model in the application case study on implementation;
Fig. 4 (b) is the Temperature Distribution cloud atlas in section at 50mm before valve.
Fig. 5 (a) is position view of the temperature measuring point before valve in the application case study on implementation;
Fig. 5 (b) is position schematic diagram on temperature measuring point section in the application case study on implementation.
Fig. 6 is leakage flow calculation flow chart in the application case study on implementation.
In figure, 1 is check-valves, and 2 be peace note pipe, and 3 be ozzle, and 4 be main pipeline
Specific embodiment
For convenience those skilled in the art understand that technology contents of the invention, illustrate the present invention in conjunction with lower example.
A kind of monitoring method of the internal drainage of check valve leakage current amount calculated based on wind-structure interaction of the present invention, including following steps
It is rapid:
The first step, the geometrical model for establishing pipeline to be analyzed, and mesh generation is carried out to fluid domain and solid domain, with meter
Fluid operator machine software carries out fluid and structural simulation, obtains valve preceding pipeline temperature field:
Channel temp t is responsible for before known pipeline position to be analyzed valvehot, main pipeline flow Q before valvehotAnd pipe internal leakage stream after valve
The temperature t of bodycoldOpereating specificationDifferent leakage flow QLeak
Tube wall temperature data before lower valve, (one group of tube wall temperature here is one to the corresponding one group of tube wall temperature of each leakage flow
Temperature data on the upper tube wall circumference of section);
Based on fluid and structural simulation, monitoring cross section temperature measuring point before n valve in the circumferential direction, the 1st temperature measuring point are determined before valve
Set on the top side of sectional position, the n-th temperature measuring point is located at the least significant end of sectional position;The temperature of two neighboring temperature measuring point successively takes
Average value is denoted as Ti, i=1,2,3 ..., n-1, the area of monitoring cross section is denoted as S between two adjacent temperature measuring pointsi, calculated according to formula (8)
Obtain the equivalent face mean temperature T of tested pointA, maximum temperature T in whole cross section is calculated according to formula (9)topWith minimum temperature
TBottomTemperature difference Δ T, the equivalent face mean temperature T of tested pointAWith temperature difference Δ T as tube wall temperature characteristic parameter,
Δ T=Ttop-TBottom (9)
In formula:
The determination of fluid and structural simulation and leakage flow correlation under second step, different flox conditions:
To fluid and structural simulation is carried out under different flox conditions, tube wall temperature feature before valve is joined using least square method
Number data are analyzed and are handled, and tube wall temperature characteristic parameter is with leakage flow Q before obtaining valveLeak, leak fluid temperature
tcold, supervisor channel temp thotAnd main pipeline flow QhotChanging rule;Tube wall temperature characteristic parameter is (to be measured before fitting valve
The equivalent face mean temperature T of pointA, temperature difference Δ T) with leakage flow QLeakQuantitative relation formula be T=f (thot,Qhot,tcold,
QLeak), T represents T in formulaA、ΔT。
According to T obtained above and QLeakMathematical relationship, solve leakage flow QLeakJoin with tube wall temperature feature before valve
Number, supervisor channel temp thotAnd main pipeline flow Qhot, leak fluid temperature tcoldBetween relational expression be QLeak=g (thot,Qhot,
tcold, T),
QLeak=g (thot,Qhot,tcold, T) and it is tube wall temperature characteristic parameter, supervisor channel temp t before valvehotAnd supervisor
Road flow Qhot, leak fluid temperature tcoldWith leakage flow QLeakBetween complicated valve leakage quantitative prediction model.
Third step carries out real-time monitoring according to prediction model:
Temperature measuring point installs thermocouple before valve, and temperature at temperature measuring point before valve is sent to monitoring server end in real time, is obtained
Thermally stratified layer temperature profile parameter, calculate leakage flow further according to valve leakage quantitative prediction model, compare by two temperature
Both the correspondence leakage flow that degree characteristic parameter obtains, if the leakage flow relative error of the two is less than leak threshold, take
Arithmetic average, if being not less than leak threshold, taking the two the larger value is final leakage stream magnitude, realizes real-time monitoring.
The method of the present invention principle is:
For the pipeline monitored, fluid and structural simulation model is established, is parameterized according to actual condition parameter area
Simulation calculates, and determines and monitors thermally stratified layer temperature in position under different supervisor's channel temps and flow, the temperature of leak fluid and leakage flow
Characteristic parameter is spent, using least square method determine thermally stratified layer temperature profile parameter and is responsible for channel temp and flow, leak fluid
Relational expression between temperature and leakage flow.When practical application, according to power plant's data (supervisor's channel temp, main pipeline flow, leakage
The temperature of fluid) and the thermally stratified layer temperature profile parameter that position monitoring obtains is monitored, it is leaked using identified relational expression
Flow.
The foundation of model, the subdivision method of grid are common sense in the field in the present invention.
Channel temp t is responsible for before the pipeline position to be analyzed valve mentioned in the first stephot, flow QhotAnd pipe internal leakage stream after valve
The temperature t of bodycoldOpereating specificationDifferent leakage flow QLeak
Tube wall temperature data are given data before lower valve, can be provided by power plant.
The pipeline that this patent method is suitable for having special thermal stratification, for the inspection of core pipe line check valve internal leakage
It surveys, the direction leaked here: to internal leakage, is particularly suitable for core pipeline, if other pipelines have similar phenomenon, this also can be used
Method is monitored.
Embodiment 1
Illustrate embodiment so that PWR main pipe line hot arc pacifies note pipeline internal drainage of check valve leakage as an example.In the upper pump that fills to high pressure
In the case that the RIS isolating valve that peace is injected back into road does not leak, there is vortex and free convection phenomenon in 1 downstream pipeline section of RCP check-valves,
Concrete condition may refer to the schematic diagram of Fig. 1.The development length of this vortex be by the pressure of primary Ioops, temperature, flow and
What the space layout of RCP check valve downstream pipeline section was determined.Pacifying note pipeline due to RCP check-valves valve downstream is high-pressure side, to the greatest extent
During pipe normal operating the valve be close, but filled on RCV pump high pressure low temperature water will enter primary Ioops by RCP check-valves,
The injection of cold water forms thermally stratified layer, and the hot fluid reciprocation penetrated into main pipeline turbulent flow in a certain range, in certain model
The interior hot and cold alternation for forming fluid is enclosed, and eventually leads to heat fatigue.So needing to monitor at the valve location whether have leakage
Occur.
By the calculating to channel temp, flow and branch pipe leak fluid temperature is responsible under different operating conditions, different leakages are acquired
Outside wall temperature under flow establishes valve leakage tube wall temperature characteristic parameter, is responsible for channel temp and flow, branch pipe (after valve
Pipe) internal leakage fluid temperature and leakage flow between complicated valve leakage quantitative prediction model.
Main pipeline flow, main pipeline fluid temperature (F.T.), leak fluid temperature and leakage flow etc. flow first under different operating conditions
Fluid-solid conjugated heat transfer calculating is carried out under the conditions of dynamic, determines the temperature of temperature tested point before valve, and the temperature for then arranging tested point is special
Parameter is levied, and is associated with flox condition, and then obtains and temperature profile parameter and the corresponding leakage stream of flox condition can be monitored
Measure correlation.
It is calculated under 1 a certain major loop flow and leakage flow based on fluid-solid conjugated heat transfer
1.1 model foundations and grid dividing
Primary Ioops main pipeline 4 hot arc bore Φ 736.6mm, wall thickness 67mm, length 4904mm are horizontally arranged.Peace note pipe 2
Outer diameter Φ 168.3mm, SCH160 (internal diameter Φ 131.78mm), are connected horizontal by 30 ° with main pipe hot leg by ozzle 3,
Check-valves 1 is installed thereon.SolidWorks establishes the geometrical model for having main pipe hot leg and peace note pipe, and the moon is established in valve
Valve leak is simulated in the aperture of thread form, as shown in Fig. 2, arrow direction is fluid flow direction in Fig. 2.
The grid dividing of fluid domain and solid domain is carried out using ANSYS ICEM software, wherein fluid domain grid and solid domain
Grid imports Fluent and carries out flow field, Simulating Calculation of Temperature Field.Fluid domain and solid domain all use six faces of higher mesh quality
Body structured grid, lattice number are 740,000, meet the requirement of grid independence.Grid is as shown in Figure 3.(a) is integral net in Fig. 3
Lattice;(b) grid of branch pipe cross section;(c) grid of main pipeline and ozzle junction;(d) grid at ozzle.
1.2 physical parameter
The material of supervisor and peace note pipeline section is Z2CN 18-10 (austenitic stainless steel).Metal material attribute is derived from RCC-M
Specification.The physical parameter of water is set as varying with temperature, will using the application program based on the physical property data in technical data
Correlation is calculated as the function of temperature.
1.3 the determination of boundary condition
The outer wall of main pipeline and Re An note pipe is set as adiabatic boundary.Ignore pipe end cross section heat output, therefore also sets
For adiabatic boundary.The interface of fluid domain and solid domain is set as coupling wall boundary condition.
Heat peace note pipe import is set as speed inlet boundary;The flow when flow of main pipeline presses Power operation determines, is set as
Speed inlet boundary condition, pressure export, pressure 15.5MPa.In order to avoid the influence that inducer simulates mixed process, disappear
Except the interference of boundary condition stream field, while reducing calculation amount, is write by the DEFINE_PROFILE of UDF, main pipeline and heat
Peace note line inlet VELOCITY DISTRIBUTION uses fully developed VELOCITY DISTRIBUTION.It will be developed completely using the secondary development UDF of fluent
VELOCITY DISTRIBUTION afterwards is added in main, branch pipe inlet boundary condition.
1.4 mathematical calculation models are established
Use Pressure-Based solver;Consider that gravity influences.Energy equation is opened, temperature field is solved: being used
Transition SST shear stress turbulence model uses StandardWall Function near wall processing mode.energy
Residual error is controlled 10-6Hereinafter, other residual error item controls are 10-4Below, it is believed that stable state calculating has had reached stable convergence
State.Using experimental data published in the document (thermal stratification in Guo Depeng, Lu Daogang, Feng Yuheng nuclear power plant pipeline
[J] atomic energy science and technology, 2007 (05): 570-574.) calculated result to be verified, the two is very identical.
The analysis of 1.5 calculated results
Along the circumferential direction equal radians are divided into temperature measuring point before n valve to monitoring cross section before valve, and the 1st temperature measuring point is set to sectional position
Top side, the n-th temperature measuring point is located at the least significant end of sectional position;The temperature of two neighboring temperature measuring point, which is successively averaged, to be denoted as
Ti, i=1,2,3 ..., n-1, the area of monitoring cross section is denoted as S between two adjacent temperature measuring pointsi, tested point is calculated according to formula (8)
Equivalent face mean temperature TA, maximum temperature T in whole cross section is calculated according to formula (9)topWith minimum temperature TBottomTemperature
Poor Δ T, the equivalent face mean temperature T of tested pointAWith temperature difference Δ T as tube wall temperature characteristic parameter,
Δ T=Ttop-TBottom (9)
In formula:
The determination of fluid and structural simulation and leakage flow correlation under 2 different flox conditions
To carrying out fluid and structural simulation under different flox conditions, using monofactor quantity method analysis primary Ioops main pipeline with
Fluid properties (the main pipeline fluid temperature (F.T.) t of peace note pipehot, main pipeline fluid flow Qhot, be connected with main pipeline peace note pipe non-return
Leak fluid temperature t after valve valve in branched pipecoldWith the leakage flow Q for passing through check-valvesLeak) to pipe outer wall temperature before check-valves
The influence of degree obtains the mathematical relationship between unitary variant and tube wall temperature based on least square method, obtains tube wall temperature
Functional relation T=f (t between characteristic parameter and influence factorhot,Qhot,tcold,QLeak)。
In formula: T represents TA, Δ T, thot、tcold、Qhot、QLeakFor regression variable (independent variable), ε is to calculate tube wall temperature spy
The random error of value indicative, α0······α5, β0······β7For required coefficient.
According to obtained T and QLeakMathematical relationship, solve leakage flow QLeakWith tube wall temperature characteristic parameter before valve,
It is responsible for channel temp thotAnd main pipeline flow Qhot, leak fluid temperature tcoldBetween relationship, so as to establish internal leakage flow QLeak
Valve leakage quantitative prediction model:
QLeak=g (thot,Qhot,tcold,T) (12)
3. carrying out real-time monitoring according to prediction model
In thot=326 DEG C, Qhot=100%QVolume、tcold=54 DEG C, QLeakUnder=400L/h operating condition, pass through fluid solid heat couple
Steady-state simulation calculates, and discovery produces apparent thermal stratification, temperature profile such as Fig. 4 (a)-Fig. 4 (b) in check-valves front end
It is shown.As shown in Fig. 5 (b), the temperature measuring point is located at the peace note pipe branched pipe of check-valves front end for temperature monitoring point position and distribution
Position on road before the check-valves to be detected at 50mm is equipped with 7 thermocouples in measurement position altogether, along the circumferential direction between
It is arranged every 30 °.Under above-mentioned operating condition, t0 °=557.282K, t30°=557.241K, t60°=553.442K, t90°=
537.567K, t120°=509.886K, t150°=468.64K, t180°=452.42K, tube wall temperature characteristic parameter determines before valve
It is as follows:
Δ T=T1-T7 (9)
In formula: SiThe area of monitoring cross section, T between two adjacent temperature measuring pointsAFor the equivalent face mean temperature of tested point, Δ T
For the temperature difference to temperature measuring point 1 and 7.
According to pipe outer wall temperature measuring point temperature, successively it is calculated by formula (1)~(6): T1=557.262K, T2=
555.342K, T3=545.505K, T4=523.727K, T5=489.263K, T6=460.53K;S1=6.3958 × 10-4m2,
S2=3.7058 × 10-3m2, S3=6.7721 × 10-3m2, S4=S3, S5=S2, S6=S1, obtained according to formula (7): S=2.2235
×10-2m2;Tube wall temperature characteristic parameter data before valve: T are calculated according to formula (8), (9)A=529.031K, Δ T=
104.862K。
The equivalent face mean temperature T of temperature measuring point is calculated according to formula (8)A is surveyed, by TA is surveyedIt is updated in formula (12), obtains QTA,
The temperature difference Δ T of the temperature measuring point is calculated according to formula (9)It surveys, by Δ TIt surveysIt is updated in formula (12), obtains QΔT, by letting out for the two
Leakage current amount relative error is compared with leak threshold, if meeting formula (13), then taking the arithmetic average formula of both
(14), if not satisfied, then taking the larger value formula (15) between two.Leakage flow calculation flow chart is as shown in Figure 6.The present embodiment
Middle leak threshold is taken as 5%.
QLeak=max { QTA,QΔT} (15)
The analysis of 4 implementation results
Multiple decision coefficients R2Reflect the fitting degree of multiple regression curve, value range is between [0,1], and R2More connect
It is bordering on 1, illustrates that regression equation is fitted better, R2Closer to 0, illustrate that regression equation is fitted poorer.Formula (10), (11)
Coefficient of determination is respectively as follows: R2=0.9968, R2=0.9957, it is in close proximity to 1, illustrates that the fitting degree of regression straight line is very good.
It is the ratio of the total sum of squares of deviations of regression sum of square Zhan.
The significance test of regression equation is to examine to whether there is one between dependent variable and all independents variable significantly
Linear relationship, also referred to as overall significance test.The method of inspection is to examine using F whether to analyze difference between the two
Significantly.Even F >=Fα(p, n-p-1) (α is significance, is taken as 0.05;P is variable number;N is the number for seeing side value),
Then there are linear relationships between dependent variable and each independent variable, if F < Fα(p, n-p-1), then between dependent variable and each independent variable not
There are linear relationships.Herein: F=1895.8275, F > F0.05(5, n-p-1)=2.21, therefore regression equation is significant, it was demonstrated that TAWith
There are linear relationships between all independents variable.F=956.2895, F > F0.05(7, n-p-1)=2.01, therefore regression equation is aobvious
It writes, it was demonstrated that there are linear relationships between Δ T and all independents variable.
By examining above, show to return obtained TA=f (thot,Qhot,tcold,QLeak), Δ T=y (thot,Qhot,
tcold,QLeak) expression formula have believable precision.Correlation result is fitted generally to coincide with numerical simulation result preferably, it is maximum
Error is 7.69% and 3.52%, and minimal error is 0.043% and 0.0438%, this illustrates correlation precision with this condition
Preferably, it also indicates that simultaneously, which is applied to have preferable reliability when the quantitative calculating of leakage flow.Wherein, error
Calculating by formula (15) calculate:
Those of ordinary skill in the art will understand that the embodiments described herein, which is to help reader, understands this hair
Bright principle, it should be understood that protection scope of the present invention is not limited to such specific embodiments and embodiments.For ability
For the technical staff in domain, the present invention can be there are many approximating method, formula form and geometrical model, grid dividing construction method
Deng.All within the spirits and principles of the present invention, any modification, equivalent replacement, improvement and so on should be included in the present invention
Scope of the claims within.
The present invention does not address place and is suitable for the prior art.
Claims (3)
1. a kind of monitoring method of the internal drainage of check valve leakage current amount calculated based on wind-structure interaction, is included the following steps:
The first step, the geometrical model for establishing pipeline to be analyzed, and mesh generation is carried out to fluid domain and solid domain, it is flowed with calculating
Mechanics software carries out fluid and structural simulation, obtains valve preceding pipeline temperature field:
Channel temp t is responsible for before known pipeline position to be analyzed valvehot, main pipeline flow Q before valvehotAnd pipe internal leakage fluid after valve
Temperature tcoldOpereating specificationDifferent leakage flow QLeakLower valve
Preceding tube wall temperature data, the corresponding one group of tube wall temperature of each leakage flow;
Based on fluid and structural simulation, determine before valve that temperature measuring point before n valve, the 1st temperature measuring point are set to monitoring cross section in the circumferential direction
The top side of sectional position, the n-th temperature measuring point are located at the least significant end of sectional position;The temperature of two neighboring temperature measuring point is successively averaged
Value is denoted as Ti, i=1,2,3 ..., n-1, the area of monitoring cross section is denoted as S between two adjacent temperature measuring pointsi, it is calculated according to formula (8)
The equivalent face mean temperature T of tested pointA, maximum temperature T in whole cross section is calculated according to formula (9)topWith minimum temperature TBottom
Temperature difference Δ T, the equivalent face mean temperature T of tested pointAWith temperature difference Δ T as tube wall temperature characteristic parameter,
Δ T=Ttop-TBottom (9)
In formula:
The determination of fluid and structural simulation and leakage flow correlation under second step, different flox conditions:
To fluid and structural simulation is carried out under different flox conditions, using least square method to tube wall temperature characteristic parameter number before valve
According to being analyzed and being handled, tube wall temperature characteristic parameter is with leakage flow Q before obtaining valveLeak, leak fluid temperature tcold, it is main
Pipe temperature thotAnd main pipeline flow QhotChanging rule;Fit tube wall temperature characteristic parameter and leakage flow before valve
QLeakQuantitative relation formula be T=f (thot,Qhot,tcold,QLeak), T represents T in formulaA,ΔT;
According to T obtained above and QLeakMathematical relationship, solve leakage flow QLeakWith tube wall temperature characteristic parameter before valve,
It is responsible for channel temp thotAnd main pipeline flow Qhot, leak fluid temperature tcoldBetween relational expression be QLeak=g (thot,Qhot,tcold,
T),
QLeak=g (thot,Qhot,tcold, T) and it is tube wall temperature characteristic parameter, supervisor channel temp t before valvehotAnd main pipeline stream
Measure Qhot, leak fluid temperature tcoldWith leakage flow QLeakBetween complicated valve leakage quantitative prediction model;
Third step carries out real-time monitoring according to prediction model:
Temperature measuring point installs thermocouple before valve, temperature at temperature measuring point before valve is sent to monitoring server end in real time, obtained heat
Stratification temperature characteristic parameter calculates leakage flow further according to valve leakage quantitative prediction model, compares by two temperature spies
The correspondence leakage flow that sign parameter obtains takes counting for the two if the leakage flow relative error of the two is less than leak threshold
Average value, if being not less than leak threshold, taking the two the larger value is final leakage stream magnitude, realizes real-time monitoring.
2. the monitoring method of the internal drainage of check valve leakage current amount according to claim 1 calculated based on wind-structure interaction,
It is characterized in that, in second step, temperature measuring point is located at before check-valves at 50mm before n valve.
3. the monitoring method of the internal drainage of check valve leakage current amount according to claim 1 calculated based on wind-structure interaction,
It is characterized in that, this method is applied to the monitoring of core pipeline peace note pipeline internal drainage of check valve leakage current amount.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910089743.3A CN109783972B (en) | 2019-01-30 | 2019-01-30 | Method for monitoring leakage flow in check valve based on fluid-solid coupling analysis and calculation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910089743.3A CN109783972B (en) | 2019-01-30 | 2019-01-30 | Method for monitoring leakage flow in check valve based on fluid-solid coupling analysis and calculation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109783972A true CN109783972A (en) | 2019-05-21 |
CN109783972B CN109783972B (en) | 2022-11-08 |
Family
ID=66503642
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910089743.3A Active CN109783972B (en) | 2019-01-30 | 2019-01-30 | Method for monitoring leakage flow in check valve based on fluid-solid coupling analysis and calculation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109783972B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110428915A (en) * | 2019-08-05 | 2019-11-08 | 辽宁红沿河核电有限公司 | A kind of nuclear power unit safety injection system minimum discharge test specific purpose tool |
CN111308337A (en) * | 2020-03-31 | 2020-06-19 | 张铭源 | Performance evaluation method for centrifugal fan inlet guide vane regulating valve |
CN112924115A (en) * | 2021-03-16 | 2021-06-08 | 中电华创(苏州)电力技术研究有限公司 | Device and method for monitoring internal leakage of high-temperature and high-pressure pipeline valve |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5637789A (en) * | 1995-09-29 | 1997-06-10 | Lawson; William J. | Fluid leak detector |
US20140283582A1 (en) * | 2011-09-27 | 2014-09-25 | Delaware Capital Formation, Inc. | Method And Apparatus For Determining The Thermal Status Of Fuel In A Line Leak Detection System |
CN104654024A (en) * | 2015-02-12 | 2015-05-27 | 常州大学 | Method for locating and analyzing leakage of city gas pipeline based on GRNN (Generalized Regression Neural Network) |
CN205090197U (en) * | 2015-07-09 | 2016-03-16 | 南京声宏毅霆网络科技有限公司 | Leak detection system and monitoring facilities of pipeline |
CN106015951A (en) * | 2016-07-07 | 2016-10-12 | 北京信息科技大学 | Gas transmission pipeline leakage detection system and method suitable for multiple state changes |
CN106650120A (en) * | 2016-12-27 | 2017-05-10 | 中国石油大学(北京) | Measuring and calculating method for leakage rate of liquid hydrocarbon pipeline |
CN108758353A (en) * | 2018-06-12 | 2018-11-06 | 青岛汇安谷科技发展有限公司 | Distribution type fiber-optic many reference amounts pipe leakage positioning alarm system and leakage locating method |
-
2019
- 2019-01-30 CN CN201910089743.3A patent/CN109783972B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5637789A (en) * | 1995-09-29 | 1997-06-10 | Lawson; William J. | Fluid leak detector |
US20140283582A1 (en) * | 2011-09-27 | 2014-09-25 | Delaware Capital Formation, Inc. | Method And Apparatus For Determining The Thermal Status Of Fuel In A Line Leak Detection System |
CN104654024A (en) * | 2015-02-12 | 2015-05-27 | 常州大学 | Method for locating and analyzing leakage of city gas pipeline based on GRNN (Generalized Regression Neural Network) |
CN205090197U (en) * | 2015-07-09 | 2016-03-16 | 南京声宏毅霆网络科技有限公司 | Leak detection system and monitoring facilities of pipeline |
CN106015951A (en) * | 2016-07-07 | 2016-10-12 | 北京信息科技大学 | Gas transmission pipeline leakage detection system and method suitable for multiple state changes |
CN106650120A (en) * | 2016-12-27 | 2017-05-10 | 中国石油大学(北京) | Measuring and calculating method for leakage rate of liquid hydrocarbon pipeline |
CN108758353A (en) * | 2018-06-12 | 2018-11-06 | 青岛汇安谷科技发展有限公司 | Distribution type fiber-optic many reference amounts pipe leakage positioning alarm system and leakage locating method |
Non-Patent Citations (3)
Title |
---|
JONG-IL PARK: "A numerical analysis on thermal stratification phenomenon by in-leakage in a branch piping", 《JOURNAL OF MECHNICAL SCIENCE AND TECHNOLOGY》 * |
段峰波等: "基于流固耦合的轴流式止回阀开启过程仿真", 《油气储运》 * |
郭德朋等: "核电厂管线中的热分层现象", 《原子能科学技术》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110428915A (en) * | 2019-08-05 | 2019-11-08 | 辽宁红沿河核电有限公司 | A kind of nuclear power unit safety injection system minimum discharge test specific purpose tool |
CN111308337A (en) * | 2020-03-31 | 2020-06-19 | 张铭源 | Performance evaluation method for centrifugal fan inlet guide vane regulating valve |
CN112924115A (en) * | 2021-03-16 | 2021-06-08 | 中电华创(苏州)电力技术研究有限公司 | Device and method for monitoring internal leakage of high-temperature and high-pressure pipeline valve |
Also Published As
Publication number | Publication date |
---|---|
CN109783972B (en) | 2022-11-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110500511B (en) | Urban non-metal pipeline leakage positioning method | |
Xu et al. | An overview of transient fault detection techniques | |
Meniconi et al. | Experimental investigation of coupled frequency and time-domain transient test–based techniques for partial blockage detection in pipelines | |
Pan et al. | FRF-based transient wave analysis for the viscoelastic parameters identification and leak detection in water-filled plastic pipes | |
CN109783972A (en) | Monitoring method based on the internal drainage of check valve leakage current amount that wind-structure interaction calculates | |
CN110197049B (en) | Non-metal pipeline leakage positioning method based on transient inverse problem | |
CN106015951A (en) | Gas transmission pipeline leakage detection system and method suitable for multiple state changes | |
Wang et al. | Experimental study on water pipeline leak using In-Pipe acoustic signal analysis and artificial neural network prediction | |
CN206130547U (en) | Gas transmission pipeline leak testing system under multiplex condition | |
CN110375924B (en) | Power plant valve leakage quantitative evaluation system | |
Wang et al. | Leak detection in pipeline systems using hydraulic methods: A review | |
Kickhofel et al. | Turbulent penetration in T-junction branch lines with leakage flow | |
Fu et al. | Detection of two-point leakages in a pipeline based on lab investigation and numerical simulation | |
Ndalila et al. | Modeling dynamic pressure of gas pipeline with single and double leakage | |
Zhang et al. | A transient-features-based diagnostic method of multi incipient cracks in pipeline systems | |
Fu et al. | Identifying two-point leakages in parallel pipelines based on flow parameter analysis | |
CN106092351A (en) | nuclear power station pipeline inner wall temperature measuring method and device | |
Sierra-Espinosa et al. | Vibration failure in admission pipe of a steam turbine due to flow instability | |
Liu et al. | Quantitative study on internal detection signal of weak magnetic stress in oil and gas pipelines based on force-magnetic noncoaxial effect | |
Shehadeh et al. | Modelling the effect of incompressible leakage patterns on rupture area in pipeline | |
Kim | Advanced numerical and experimental transient modelling of water and gas pipeline flows incorporating distributed and local effects. | |
Yan et al. | Valve failure detection of the long-distance district heating pipeline by hydraulic oscillation recognition: A numerical approach | |
Rivas et al. | Application of mass balance and transient flow modeling for leak detection in liquid pipelines | |
Covas et al. | Water pipe system diagnosis by transient pressure signals | |
Guo et al. | Leak detection in pipelines by exclusively frequency domain method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |